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The Maxwell system in waveguides with several ends 

B. A. Plamenevskii, A. S. Poretckii Q 

A waveguide coincides with a domain G in M 3 having finitely many cylindrical outlets 
to infinity; the boundary dG is smooth. In G, we consider the stationary Maxwell system with 
spectral parameter k G M and identity matrices of dielectric and magnetic permittivity. The 
boundary dG is supposed to be perfectly conductive. In the presence of charges and currents 
we investigate the solvability of the corresponding boundary value problem supplemented 
with "intrinsic" radiation conditions at infinity. For all k in the continuous spectrum of 
the problem (including the thresholds and eigenvalues), we describe a basis in the space of 
continuous spectrum eigenfunctions, define the scattering matrix, and prove it is unitary. To 
this end, we extend the Maxwell system to an elliptic one and study the latter in detail. 
The information on the Maxwell boundary value problem comes from that obtained for the 
^ ' elliptic problem. 

1 Introduction. Formulation of the results 

Let G be a domain in M 3 , coinciding outside a large ball with the union U • • ■ U 11+ of 
finitely many non-overlapping semicylinders 

n q + = {(y«,t«):y*eW, i« > 0}, 

in ; 

CN ! where (y q ,t q ) are local coordinates in U q + and Vl q is a bounded domain in M 2 . We consider 
the Maxwell system 

O 1 i rot u 2 (x) — ku 1 (x) = / x (x), 

-idivu 2 (x) = h\x), (1.1) 

— i rot u (x) — ku 2 (x) = f 2 (x), 

^| idivu 1 (x) = h 2 (x) 
H ■ 

in G with boundary conditions 

u(x) x u\x) = 0, (u 2 (x), u(x)) = 0, x E dG; (1.2) 

here u l and u 2 are vector valued functions with three components (the electric and magnetic 
vectors respectively), k is a spectral parameter, v is the outward normal to dG, (u,u) are 
i/XM are scalar and vector product, and p , h? are given functions. The boundary conditions 
correspond to perfectly conductive boundary, which means that at dG there vanish the 
tangent component of u 1 and the normal one of u 2 . The system (II. ID is overdetermined 
(eight equations and only six indeterminate functions). The compatibility conditions 

div/ x (x) - ikh 2 (x) = 0, x e G, 

divf 2 (x) + ikh l (x) = 0, x EG, (1.3) 
{f 2 (x),v(x)) = 0, x G dG 
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are necessary for the solvability of problem (jl.ip . ( II. 2p . 

Preparatory to describing the continuous spectrum eigenfunctions and the radiation 
conditions, we first associate with problem ( II. ip . (II. 2p some operator pencils and then define 
" waves " needed for asymptotic formulas. Let us write down (11.11) . (II. 2p in the form 

M{D,k)U(x) = JF(x), xeG; u 1 (x) T = 0, u 2 (x) u = 0, x G dG, (1.4) 

where D = (D\, D 2 , D3), Dj = —id/dxj, U = (u 1 ,^ 2 ), while u 1 (x) T and u 2 (x) v denote the 
tangent and normal component of u l (x) and u 2 {x). 

Let us consider problem (jl.4p in the cylinder G := Vt x R = {x = (xi,x 2 ,x 3 ) : 
(xi,x 2 ) G ^,^3 G M}, where f2 is a bounded domain in IR 2 with smooth boundary dQ. 
For the vectors $ = (<p,ip) with components ^> in C°°(f2;C 3 ) satisfying the boundary 
conditions (p T = 0, ip„ = at <9f2 (we assume here that the outward normal to dQ is the 
vector v — (ui, u 2 , 0)), we define an operator pencil C 3 A H- 3H(A, k) by the equality 

9Jt(A, k)$(xi,x 2 ) = exp (-iAx 3 )A4(.D, fc)(exp (i\x 3 )$(x 1 ,x 2 )). (1.5) 

A number A G C is said to be an eigenvalue of the pencil 9Jt(-, fc) if there exists a vector $ 7^ 
such that 2Jt(A, fc)<& = 0; then $ is called an eigenvector of the pencil. The eigenvalues of 
9Jt(-, fc) are resting on the coordinate axes of the complex plane, accumulate only at infinity, 
while the real axis contains at most finitely many eigenvalues. We assume that the domain 
Q is 1-connected. A value k 7^ is called a threshold if A = turns out to be an eigenvalue 
for Wl(-,k). The set of thresholds is infinite, symmetric about the coordinate origin, and 
accumulates only at infinity. The first threshold on the semiaxis < k < 00 is equal to JJIq, 
where /x is the minimal positive eigenvalue of the Neumann problem 

Au(x) + fiu(x) =0, if!!, d u u(x) = 0, x G dQ, 

A being the Laplace operator in Q. Let k' and k" are neighboring thresholds and let < 
k' < k" . The sum of multiplicities n{k) of the real eigenvalues A of the pencil 97t(-, k) is even 
and constant for \k\ G (k',k"); to every eigenvalue A there correspond only eigenvectors, 
while generalized eigenvectors do not exist. We denote by Ao and $0 a rea l eigenvalue and a 
corresponding eigenvector of k). The function 

G 3 (y, t) ^ P{y, t) = exp {iX t)^ {y) (1.6) 

satisfies the homogeneous problem (II. 4p in the cylinder G — {x — (y, t) : y G fl, t G R}. Each 
of the solutions is a linear combination of TE- and TM-modes (see, e.g., [T]). The dimension 
of the linear space of solutions of the form ( \1.6h is equal to n(k). One can choose a basis in 
the space subject to certain orthogonality and normalization conditions so that half of the 
basis consists of "waves outgoing to +00" , and the other half consists of "waves incoming 
from +00" . 

Let us turn to the "general" domain G with several cylindrical ends . . . , LT^. We 
introduce the pencil 0Jt <? (- , k) in II 9 = Q q x K., q = 1, . . . , N, denote by X 9 the set of thresholds 
for Wl q , and define X := X 1 U ■ • • U % N . Assume now that k', k" G X, < kl < k", while 

(k', k") n X = 0. We also suppose that x(k) := H X {k) H h x N (k), where k G (k 1 , k") and 

x q {k) is the sum of multiplicities of the real eigenvalues of the pencil DJl q (-,k). The value 
x(k) is constant for |fc| G (k', k") and equal to an even number 2T. 
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Denote by \ a function in C°°(R) such that < x(t) < 1? x{t) = 1 for t > T, and 
x(t) =0 for t < T — 1, where T is sufficiently large number. In local coordinates of 11+ , 
introduce the function (y q ,t q ) i— > x(t 9 )P(i/' ? , t 9 ), P being a solution of the form (jl.6p to the 
homogeneous problem (jl.4p in II 9 . We can assume that suppx-P C G D II + for a large T 
in the definition of x- Extending %P by zero to G, we obtain a smooth function given in 
G. Taking for P in II 9 a wave incoming from +oo (outgoing to +oo), we obtain a function 
called incoming (outgoing) wave in the domain G. In such a way we define incoming waves 
wf, . . . , and outgoing waves w^, . . . , w^- in G. 

As before, assume that \k\ G (&', fc"). A smooth bounded in G vector- valued function 
U = (w 1 ,?/ 2 ) is called a continuous spectrum eigenfunction (CSE) of problem (II. ip . (jl.2p if 
W satisfies the homogeneous problem ( II. ip . ( II. 2p and does not belong to the space L 2 (G). 

Denote by p$ a smooth positive function on G, given on U q + nG by p$(y q , t q ) = exp (5t q ) 
with q = 1, . . . , N and a positive 5. We choose such a 5 sufficiently small so that the strip 
{A G C : | Im A | < 5} contains no eigenvalues of the pencils DJl q (-, k) except for the real ones. 
Introduce the space Wg(G) of functions in G with norm 

\\u;W>(G)\\ :=(E / \D°{psu)\ 3 dz) 1 * 

for I = 0, 1, ... . A number k is called an eigenvalue for problem ( II. ip . (II. 2p if there exists 
in L 2 (G) a (smooth) solution to the homogeneous problem; as a rule, we do not distinct 
in notations the space of scalar and vector valued functions. Every eigenvalue is real; the 
dimension of eigenspace is always finite. The eigenvalues can accumulate only at infinity. 
Any eigenfunction turns out to be in W\{G) for each I. 

Denote by kerA^ (P, k) the eigenspace of problem (ll.ip . ( jl.2[) (possibly, trivial) corresponding 
to a number k. For the sake of formulation simplicity we suppose in Theorems 11.11 and 11.21 
as before, that k', k" G X, < kl < k", and (k', k") n % = 0. 

Theorem 1.1. Let \k\ G (k', k") and let 2T be the sum of multiplicities of all real eigenvalues 
of the pencils 07t 9 (• , k), q — 1, . . . , N. Then there exist solutions W/(-, k), j = 1, . . . , T, to 
the homogeneous problem ( II. ip . ( II .2p such that 

x 

k) - w+(; k)-"£ s jq {k)w-{; k) g W l 5 (G), 

q=l 

where / = 1, 2, . . . If k is not an eigenvalue of problem (II. ip . ( II. 2p . then the solution W/( - , k) 
with above property is unique and (■, k), . . . , Wy(-, k) form a basis in the space of CSE 
corresponding to k. If k turns out to be an eigenvalue of problem ( II. ip . ( II. 2p . then Wj(-,k) 
is determined up to an arbitrary term in kerAi(D, k) and any (•, k), . . . , (■, k) form 
a basis modulo kerA^(P, k) in the space of CSE. 
The matrix 

s (k) = (s jq (k))l q=1 

is unitary. It is independent of the choice ofWj~(-,k) in the case that k is an eigenvalue of 
problem (fTTp . (Oj) . 
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The theorem can be generalized for the thresholds as well. The matrix s is called the 
scattering matrix. It is defined for k satisfying k 2 > /%t, where /igtjj = minj/iQ, . . . , /iff } and 
(/ip) 1 / 2 is the first positive threshold for the pencil OTt- 7 in fP. In fact, there is no wave for 
k 2 < H<m that could transfer energy. 

The incoming and outgoing waves can "exchange roles" . In particular, there exist 
solutions W~(-, k) to the homogeneous problem (II. ip . f ll.2jl such that 

T 

k) - wj(; k) - 5>*(*K(-, k) e W l s {G). 

3=1 

The matrices (sj q (k))J q=1 and (tj q (k))J q=1 are mutually inverse. 

The solvability of the problem (11.1)) . (11.21) supplemented with intrinsic radiation 
conditions is established by the following theorem. 

Theorem 1.2. Let \k\ G (k',k") and let Ci>--->Cm be a basis in the space kerM(D, k) 
of eigenvectors of problem (jl.ip . (1 1 . 2 p . VKe a/so assume that J 7 = (f 1 , h 1 , f 2 , h 2 ) is in 
Wg~ 1 (G, C 8 ), (I > 1), satisfies the compatidility conditions (II. 3p and moreover (/, Q)q = 
for j = 1, . . . , m, where f = (f 1 , f 2 ) and (•, -)q is the inner product on L 2 (G). Then there 
exists a solutionis = (m 1 ,^ 2 ) with the radition conditions 

V:=U- ci^r - . . . c T w T G W l 5 (G; C 6 ), 

where Cj = i(f, Y\?J)g- Such a solution U is determined up to an arbitrary term in ker M. (D, k) 
and there holds the inequality 

||V;^(G;C 6 )|| + | Cl | + ■ •• + |c T | < constdlJ 7 ; W l 5 ~ 1 (G, C 8 )|| + \\p s V; L 2 (G, C 6 ) ||). (1.7) 

A solution Uq that satisfies the additional conditions (Uq, Q)g — is unique and there holds 
the estimate (II. 7p with right-hand changed for const|| J 7 ; W l s ~ l {G, C 8 ) || . 

The study of problem (II. ip . (II. 2p begins with extension of the overdetermined Maxwell 
system to an elliptic system. To this end we use the orthogonal extension method suggested 
by I.S. Gudovich, S.G. Krein, and I.M. Kulikov (see e.g. [2], [3] and references there). As a 
result, there arises an elliptic boundary value problem self-adjoint with respect to a Green 
formula. The general problems of this type in domains with cylindrical outlets to infinity were 
studied in j3]. In particular, the intrinsic radiation conditions were described, the solvability 
of the boundary value problem with those radiation conditions was established, the unitary 
scattering matrix was introduced. When analyzing the obtained elliptic problem, we clarify 
its specific properties coming from the Maxwell system. To this end we investigate in detail 
the operator pencil generated by the elliptic problem. Then we derive the information on 
the Maxwell system from that obtained for the elliptic one. 

From numerous mathematical works devoted to the Maxwell system in waveguides 
we set off two lines of investigation. One of the lines is related to the Wiener- Hopf technique 
and the mode matching method. Surveys of the methods are given in [5], [B]. The other line 
is presented in [7j (see also references therein). The methodology in these works is connected 
with cylindrical waveguides and dielectric and magnetic permittivity independent of the axial 
variable. 
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We use neither the methods nor the results of the works mentioned in the preceding 
paragraph. The elliptic extension of Maxwell system provides all the advantages of elliptic 
situation, in particular, the possibility of localization, freedom in choosing waveguide geometry. 
In this paper, we consider waveguides with identity matrices of dielectric and magnetic 
permittivity; in another paper, we are going to show that our approach also suggests reasonably 
wide freedom in choosing waveguide medium. 

2 Augmented Maxwell system 

2.1 Elliptic boundary value problem 

We now pass on to the "orthogonal extension" of system (II. ip ([2], see also [3]). Namely, in 
the domain G we introduce the boundary value problem 

i rot u 2 (x) + i Va 2 (i) — ku 1 ^) = 

-idivu 2 (x) - ka\x) = h\x), (2.1) 
— i rot n 1 (a;) — i Va 1 (x) — ku 2 (x) = f 2 (x), 
i div u (x) — ka 2 (x) = h 2 (x) 

with boundary conditions 

u(x) x u l (x) = g 1 (x), (u 2 (x),v(x)) = g 2 {x), a 2 (x) = g 3 (x), x G dG; (2.2) 

here u l ,u 2 are vector valued functions with three components and a 1 , a 2 stand for scalar 
functions in G. Problem (12. ip . (12. 2p is elliptic. Rewrite it in the form 

A(D,k)U(x) = T(x) 1 xeG, (2.3) 
BU{x) = Q(x), x e dG, 

where D = (D 1 ,D 2 ,-D 3 ), Dj = —id/dxj, U = (u 1 , a 1 , u 2 , a 2 ). The Green formula holds 

(A(D, k)U, V) G + (BU, QV)qg = (U, A(D, k)V) G + (QU, BV) dG , (2.4) 

with U = (u 1 , a 1 , u 2 , a 2 ), V = (v 1 , b 1 , v 2 , b 2 ), and 

BU = (ux u\ (u 2 , u), a 2 ), QV = (-iv 2 , -ib\ (iv\ v)). 

The operator of problem ( 12. 3 p is self adjoint with respect to the Green formula ( 12. 4p . 

2.2 Elliptic and Maxwell operator pencils 

We consider the operator {A(D),B} of problem (12. 3p in the cylinder Q x R = {x = 
(xi,X2,Xs) : (xi,X2) G fi,X3 G M.}, where Q is a bounded domain in IR 2 with smooth 
boundary dQ. Let $ = (ip , a , ip , /3) be a vector with components tp, if; in C°°(r2;C 3 ) and a, 
(3 in C°°(0;C) satisfying 

B$ = (i/x^<V^),/3) = (2.5) 
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on dQ, where v is the outward normal to dfl; we assume that v is of the form (z/ l7 i/ 2 ,0). 
Denote by u T and u v the tangent and normal components of u on dQ and rewrite (12. 5p as 

Vr = o, i) v = o, p\an = o. 

For the vectors <3> with such properties we define the operator pencil C 3 A H- 21(A), 



2l(A)$(a;i, x 2 ) = exp (-i\x 3 )A(D)(exp (i\x 3 )®(xi,x 2 )). 



(2.6) 



For the usual operations V, rot, div, and A in Q x 
rot(A), div(A) and A(A) by 



introduce in Q the operations V(A), 



V(A)a(xi, x 2 ) = exp (— -iAx 3 )V (exp (i\x 3 )a(xi, x 2 )), 
rot(A)</9(xi, x 2 ) = exp (— zAa^rot (exp (i\x 3 )ip(xi, x 2 )), 

etc. Formulas for the usual operations can immediately be modified for the operations 
with parameter. For example, from rot rot = Vdiv — A it follows that rot(A)rot(A) = 
V(A) div(A) - A(A). For tp, ip in C°°(0; C 3 ) and a in C°°(0; C) we have 



(V(A)a,p) n 
(rot(A)^,^) n 



(a, (<p, u)) m - (a, div(X)p) n , 
(ip,ip x u) m + (<p,rot(A)^) n . 



(2.7) 
(2.8) 



Denote by H l (Q; C 8 ), / = 0, 1, . . . , the space of vectors with eight components in the Sobolev 
space H l (Q; C) of functions in Q. We write the elements $ G H l (Q; C 8 ) as $ = (<^>, a, ip, (3), 
where (p,ip E #'(fi; C 3 ) and a, /3 £ H l (n; C). For / = 1, 2, . . . set 



Dtf(fi) = {$ G #'(fi; C 8 ) : <p T = 0, ^ = 0, = 0}. 



(2.9) 



Let us consider the operator 21(A) given by (12. 6ft on the domain VH l (Q). According to the 
general theory of elliptic operator pencils (see [5]), for all A G C with the exception of 
some isolated points, the mapping 21(A) : T>H l (Q) — > H l ~ 1 (Q;C 8 )is an isomorphism. The 
mentioned isolated points are the eigenvalues of the pencil A y 21(A) of finite algebraic 
multiplicity. The components of eigenvectors and generalized eigenvectors are smooth in Q. 
For U = (if, a, tp, /3) G VH l {Vt) from fllTT]) and (gSJ) it follows that 



21(A): 



a 



/ irot(A)^ + iV(A)/3- kip \ 
—i div(A) ip — ka 
—i rot(A) p — i V(A)a — kip 
i div(A) <p — k/3 



(2.10) 



The pencil 21 is called elliptic and its restriction to {U G T>H (Q) : U = (tp, 0,^,0)} will 
be called the Maxwell pencil and denoted by Wl. The number Ao is an eigenvalue of the 
pencil SUT(- , k) if there exists a smooth nonzero vector $ = (ip, 0, ip, 0) that is subject to the 
boundary conditions p T = 0, ip v — on dfl and satisfies fc)$ = 0. 
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2.3 Eigenvalues and eigenvectors of the pencils 21 and 9JI 

Proposition 2.1. Let X be an eigenvalue of the pencil 2t(-, k) and (tp, a, ip, (3) an eigenvector 
corresponding to the eigenvalue, if = (<pi,<p2,f3) andip = (^1,^2,^3)- Assume thatk 2 — X 2 7^ 
0. Then 

A(X)a + k 2 a = in fl, d u a = on dQ, (2.11) 
A(X)(3 + k 2 (3 = OinQ, (3 = on dtt, (2.12) 
A(A)^ 3 + fcVa = inn, ^ = on d£l, (2.13) 
A(A)^ 3 + k 2 ip 3 = mVl, d u 4> z = on dtt, (2.14) 

while ifj, ipj for j = 1,2 are defined by 

(pi = (k 2 — \ 2 )~ 1 [i\di(f3 + ikd 2 ip3 — i\d 2 a + ikdi(3], 

V? 2 = (k 2 - X 2 )- 1 ^^ - ikd^s + i\d ia + ikd 2 /3], (2.15) 

^ = (k 2 - X 2 )- 1 [-ikd 2 <f3 + i\d 1 'ilj 3 -ikd 1 a-iXd 2 l3], 

il) 2 = (k 2 — X 2 )~ 1 [ikdi(ps + iXd 2 ips — ikd 2 a + iXdi(3]. 

Conversely, any nonzero vector (ip,a,ip, (3) with components satisfying 02. lip - (I2.14p and 
(I2.15P is an eigenvector of$l(-,k) corresponding to X. 

If a number X (such that k 2 ' — X 2 ^ 0) is an eigenvalue for one of the pencils 2l(- , k) and 
9Jt(-, k), then it is an eigenvalue for the other one; moreover, X turns out to be an eigenvalue 
of the pencils if and only if it is an eigenvalue for at least one of problems 02. lip and 02.12p . 
There hold the equalities 

>%(A, k) = 2xgjt(A, k) = 2x v (X, k) + 2x^(A, k), 

where x%{X, k) and >%t(A, k) are the geometric multiplicities of the eigenvalue X for the 
pencils 2l(-,/c) and ^K{-,k), while xj^(X,k) and >Cx>(X,k) are those for problems 02. lip and 

(EED. 

The case k 2 - X 2 = has been described in the next proposition. 

Proposition 2.2. Let Q be a 1-connected domain. Then: 

1. If X 2 = k 2 ^ 0, then X is an eigenvalue o/2l(-, k) and the corresponding eigenspace 
is one- dimensional and spanned by the vector $ = (ip, a, ip, (3) with components 

<p = 0, a = const 7^ 0, ipi = ip 2 = 0, ^3 = (X/k)a, (3 = 0. 

The vector $ does not belong to the domain 9Jt(- , k), and X is not an eigenvalue for 0Jt(- , k). 

2. If k = 0, then X = is an eigenvalue o/2l(-,fc) with eigenspace spanned by the 
vectors $ = (0, a, ip, (3) and $ = (0, a, ip, [3), where 

= 0, a = const 7^ 0, ip = 0, (3 = 0, 

(p = 0, a = 0, = -02 = 0, ^3 = const 7^ 0, (3 = 0. 

The vector $ does not belong to the domain of9Jl(-, k), while $ is an eigenvector for 93T(-, k). 

The generalized eigenvectors of 2t(-, k) exist if and only if is a threshold, i.e. A = 
is an eigenvalue of 2l(-, k) while k 7^ 0. We do not dwell on describing such vectors. 
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2.4 Continuous spectrum eigenfunctions. Scattering matrix 

Assume that k is not a threshold for the elliptic pencils 21 1 , . . . %L N . For every q — 1, . . . , N in 
the cylinder II 9 = Q q x R, we define solutions of the form (11. 6p to the homogeneous elliptic 
problem (12. ip . (12. 2p (to this end we can use Prop osit ions 12.11 and 12 . 2 p . The dimension of linear 
space of such solutions in IF is equal to the sum £%(k) of all real eigenvalue multiplicities 
for 2l g (-,/c); the number E%(k) is even. In the space there exists a basis satisfying proper 
orthogonality and normalization conditions so that one half of the basis consists of of "waves 
outgoing to +oo" , and the other half consists of "waves incoming from +oo" . The waves in 
the domain G 

v+, v~, j = l,...,T:=(El l (k) + ..- + I${k))/2, (2.16) 

can now be introduced like those in Section 1 after (II. 6p . 

A smooth bounded in G vector- valued function U = (u 1 , a 1 , u 2 , a 2 ) is called a continuous 
spectrum eigenfunction (CSE) of problem (12.11) . (12.21) if U satisfies the homogeneous problem 
(12. ip . (12. 2p and does not belong to the space L 2 {G). A number k is called an eigenvalue for 
problem (12. ip . (12. 2p if there exists in L 2 {G) a (smooth) solution to the homogeneous problem. 
Every eigenvalue is real; the dimension of eigenspace is always finite. The eigenvalues can 
accumulate only at infinity. Any eigenfunction turns out to be in W\(G) for each I. Denote 
by kerA(D, k) the eigenspace of problem (12. ip . ( 12. 2 p (possibly, trivial) corresponding to a 
number k. 

Theorem 2.3. Let 2T be the sum of all real eigenvalue multiplicities of the pencils 2l g (-, k), 
q = 1, . . . , N . Then there exist solutions V,- (•, k), j = 1, . . . , T, to the homogeneous problem 
( 12~TD . ( |2T2|) such that 

T 

V+(-, k) - vf(-, k)-J2 °-*,(*K(-> k) e W l s (G), 

9=1 

where I = 1, 2, . . . If k is not an eigenvalue of problem (12. ip . (12. 2p . then the solution Vj~(-, k) 
with above property is unique and Vi(-,k),...,V^(-,k) form a basis in the space of CSE 
corresponding to k. If k turns out to be an eigenvalue of problem ( 12. ip . ( 12. 2 p . then Vj~(-,k) 
is determined up to an arbitrary term in kei A(D, k) and any (■, k), . . . , V^(-, k) form a 
basis modulo kei A(D, k) in the space of CSE. 
The matrix 

a(k) = (o- jq (k))l g=1 

is unitary. It is independent of the choice o/V J + (-, k) in the case that k is an eigenvalue of 
problem (I2.ip . ( I2.2p . There exist solutions V~(-,k) to the homogeneous problem \2. \2.^ 
such that 

T 

V-(; k) - V-(; k)-J2 T jq (k)v+(; k) G W\{G). (2.17) 
9=1 

The matrices {o~j q (k))J q=1 and {Tj q (k))J q=l are mutually inverse. 

This theorem can be generalized for any real k (including the thresholds), the o~(k) 
and r{k) are defined for all k £ R, i.e., the continuous spectrum of problem (12. ip . (12.21) 
coincides with R. 
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2.5 Radiation principle 

Theorem 2.4. Let z u . . . , z d be a basis in kerA(D, k) and {J 7 , Q} e W\~ X {G) x W l f 1/2 (dG), 
while 

(J 7 , Zj) G + ((?, Qzj) dG = 0, j = l,...,d, 

where Q is the operator in the Green formula (12. 4p . Then there exists a solution U of the 
problem (12. ip . (12. 2p determined up to an arbitrary term in ker.A(.D, k) such that 

V:=U-CiV^ c t VtEW 1 s (G) : (2.18) 

where 

Cj =i(F,V-) G + i(g,QVr) dG , j = l,...,T. (2.19) 
There holds the inequality 

||V;Wj(G)|| + |c 1 | + ..- + |o r | (2.20) 

< constdl^;^ 1 ^)!! + l&W^dGW + \\p s V; L 2 (G)\\). 

A solution Uq subject to (Uq, Zj) G = 0, j = 1, . . . , d, is unique and satisfies (I2.20p with 
right-hand side replaced by const(||.F; + \\Q\ W l f 1/2 dG\\). 



3 Coming back to the non-augmented Maxwell system 

To simplify formulations, we suppose that k ^ and k is not a threshold. We will consider 
the elliptic system 

i rot u 2 (x) + i Va 2 (x) — ku (x) = f (x), 

— i div u 2 (x) — ka 1 (x) = h (x), (3-1) 
— i rot n 1 (a;) — i Va 1 (x) — ku 2 (x) = f 2 (x), 
i divw 1 (x) — ka 2 (x) = h 2 (x) 

in G with homogeneous boundary conditions 

u(x) x n 1 ^) = 0, (u 2 (x), v(x)) = 0, a 2 {x) = 0, x G <9G. (3.2) 

If J 7 = (J 1 , /i 1 , / 2 , /i 1 ) belongs to Wj _1 (G f ) with / > 2, then the compatibility conditions (11731) 
can be understood directly. In fact, such conditions can be interpreted in some generalized 
form for I = 1 as well. 

Proposition 3.1. Assume that the vector J 7 = (f 1 , h 1 , f 2 , h 2 ) belongs to W l s ~ 1 {G) and is 
subject to compatibility conditions (II. 3p . Let U = (u 1 , a 1 , u 2 , a 2 ) satisfy the problem (13. ip . 
(13. 2 p and £ne radiation condition (I2.18p . Then a 1 is a solution to the problem 

(A + k 2 )a l (x) = 0, xGG, <9„a 1 (x) = 0, x E dG, (3.3) 

and satisfies the intrinsic radiation conditions defined for the problem (13. 3p . tu/wZe a 2 zs a 
solution to the problem 

(A + k 2 )a 2 (x) = 0, xeG, a 2 (x) = 0, x E dG, (3.4) 

and satisfies the intrinsic radiation conditions defined for problem (13 ,4p . 77ms a 1 (a 2 j can 6e 
nonzero only if it is an eigenfunction of problem (13. 3p (problem (13. 4p ). 
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Let us here restrict ourselves to considering the case that k is an eigenvalue of neither 
of problems ( 13. 3p and ( 13. 4p . Then the solution U in Proposition 13. II in fact satisfies ( 13. ip and 
( 13. 2 p with a 1 = a 2 = 0. It remains to discuss the radiation conditions (I2.18p . 

The collection of waves in G needed for problem ( 13. ip . (I3.2p in Theorems 12.31 and 12.41 
consists of the two parts 

E = {ef}J =1 and r = {tf}J+?. (3.5) 

The outgoing and incoming waves are generated by eigenvectors of 9Jl q (-,k), while the 
waves 7^ are generated by eigenvectors of 2l p (-, k) which belong to none of the domains of 
pencils Wl q (-,k), q = 1, . . . , N. 

For the solution IA in Proposition 13. 1[ the radiation conditions ( 12 . 1 8[) takes the form 

U-c x e{ c r e^eW l s(G); (3.6) 

however the elliptic problem yet manifests itself in the coefficients Cj = i(J-,V^)g in (13. 6p . 
where VJ is a solution to the homogeneous problem (12. ip . (12. 2p . We have to show that the 
role of VJ is in fact played here by a solution to the homogeneous Maxwell system. It can 
be done by using the following 

Proposition 3.2. There exists a unique solution £~J = (w 1 , 0, w 2 , 0) of the homogeneous 
problem (J27TJ), (Q such that 

T 

£+(-, k) - e+(-, k)-J2 s jq (k)e-(-, k) e W}{G) (3.7) 

9=1 

for j = 1, . . . , T. The functions Wj + := (it 1 , u 2 ) consisting of the components of £^ form a 
basis in the space of continuous spectrum eigenf unctions of the Maxwell problem (jl.ip . ( II. 2p . 

In order to see that the scattering matrix a in Theorem 12.31 is block diagonal, we now 
combine Proposition 13.21 and the next 

Proposition 3.3. There exists a unique solution Qj of the homogeneous problem (I2.ip . 
(T2T2]) . such that 

T+N 

v iP% &W s(G) (3.8) 

P =i 

for j = 1, . . . , T + N. The functions £+, . . . , . . . , {?y+tv f orm a basis in the space 

of continuous spectrum eigenf unctions of problem (12. ip . ( 12. 2p . The equality a = diag(s, v) 
holds, where a = o~(k) is the scattering matrix in Theorem \2.!% 

s = s(k) = (s jq (k))l q=1 , v = v(k) = (v jq (k))J+». 

Since the matrix a is unitary, every block s and v is unitary as well. 
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